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The size of the proton - closing in on the radius puzzle 
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Abstract. We analyze the recent electron-proton scattering data from Mainz using a dispersive framework 
that respects the constraints from analyticity and unitarity on the nucleon structure. We also perform a 
continued fraction analysis of these data. We find a small electric proton charge radius, r v E = 0.841q'oi fm, 
consistent with the recent determination from muonic hydrogen measurements and earlier dispersive analy- 
ses. We also extract the proton magnetic radius, r 
based on dispersion relations. 
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0.86_ ' 03 fm, consistent with earlier determinations 
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PACS. 13.40. Gp , 14.20.Dh, 11.55.F\ 



The proton charge radius is a fundamental quantity 
of physics. It is truly remarkable that despite decade long 
experimental and theoretical efforts its precise value is not 
yet determined. The recent controversy about the size of 
the proton was triggered by the precision measurement of 
the Lamb shift in muonic hydrogen that led to a "small" 



charge radius, 



p 



0.84184(67) fm pQ. This result came 



as a big surprise as it was in stark contrast to the com- 
monly accepted "large" CODATA value of 0.8768(69) fm 
[2] , based on measurements of the Lamb shift in electronic 
hydrogen and the analysis of electron-proton scattering 
data. The large value was further strengthened by the 
high precision electron-proton scattering measurements at 
\ MAMI-C PUS]. The analysis of these data including two- 
.£h ' P no ton corrections led to r v E = 0.876(8) fm. These au- 
' thors also found a magnetic radius of the proton that 
5^ , came out much smaller than commonly accepted values, 
_C3_ r p M = 0.803(17) fm. 

On the theoretical side, a precise ab initio calculation 
based on lattice QCD is not yet available due to various 
conceptual problems to be overcome like the treatment 
of disconnected contributions that feature importantly in 
isoscalar quantities. See, e.g., [5] for a recent work. A dif- 
ferent framework, that was pioneered by Hohler and col- 
laborators a long time ago, are dispersion relations (DRs) 
for the nucleon electromagnetic form factors [6]. In such 
type of approach, smaller radii were always favored, and 
the most recent and sophisticated DR calculation led to 
values consistent with the ones from muonic hydrogen, 
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+0.008 . 



netic radius was 



004 fm [7]. The value for the proton's mag- 
0.854 ± 0.005 fm, consistent with, 
e.g., the determination based on continued fractions by 
Sick [3]. Note, however, that the conformal mapping tech- 
nique to analyze the nucleon form factors led to a large 
value of the proton charge radius [DJ. 

Given this puzzling situation, in this Letter we will 
re-analyze the MAMI data of Bernauer et al. [3] using 
dispersion relations. Our main focus will be on the cor- 
rect treatment of the analytical structure of the nucleon 
form factors that is driven by the two-pion continuum. Its 
important role was already stressed by Frazer and Fulco, 
who were able to predict the p-resonance and its influ- 
ence on the nucleons' structure a long time ago [ID) . What 
has often been overlooked since this seminal work was the 
large enhancement of the two-pion continuum on the left 
wing of the p-resonance due to a close-by pole on the sec- 
ond Riemann sheet in the elastic pion-nucleon scattering 
amplitude. This enhancement amounts for roughly half 
of the nucleon isovector size [B] . This important effect is 
also recovered in chiral perturbation theory, the effective 
field theory of QCD at low energies [IT]. We consider it 
therefore of utmost importance to include this effect in the 
re-analysis of the MAMI data. Of course, in light of the 
muon g — 2 measurement at BNL [12] , one might speculate 
about the influence of new physics as the muon data are 
more sensitive to such effects, but first one has to exclude 
possible conventional explanations - and we will offer such 
a possibility here. 
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The nucleon matrix element of the electromagnetic 
current 



{N(p')\f™\N(p)) 

= ieu(p')(j li F 1 (t) 



+ i 



2m,N 



*!»(*) «(P)« (!) 



is parameterized in terms of the Dirac, Fi(t), and Pauli, 
i*2(*)) form factors, with t = (p' — p) 2 = — Q 2 the in- 
variant momentum transfer squared, and win is the nu- 
cleon mass. For electron-nucleon scattering, we have Q 2 > 
0. Ff^ n (0) and F% "(0) are given in terms of the pro- 
ton/neutron electric charge and anomalous magnetic mo- 
ment, respectively. For the later analysis, it is useful to 
separate the form factors in their isoscalar and the isovec- 
tor parts, F? = (if + F?)/2 and F? = (F[ - Ff)/2 
for i = 1,2, correspondingly. The differential cross section 
is given most compactly in terms of the Sachs form fac- 
tors G E {t) = F 1 {t)~rF 2 (t), G M {t) = F 1 (t) + F 2 (t), with 
r = — t/4m N , so that 



da 
dQ 



da 
dQ 



GUQ 2 



(2) 



Matt e ( X + 

r) tan 2 (6*/2)] _1 is the virtual photon 



where e = [1 + 2(1 
polarization, is the electron scattering angle in the lab- 
oratory frame, and (da / dQ)Mott is the Mott cross section, 
which corresponds to scattering on a point-like particle. 
The electric and magnetic radii of the proton, that are in 
the focus of this Letter, are given by 



-6 dG E MQ 2 ) 



E,M 



g e ,m(P) d Q 2 



1/2 



(3) 



Q 2 =o, 



To analyze the cross section data from Mainz based 
on Eq. @, we use unsubtracted dispersion relations for 
the nucleon form factors. For a generic form factor in the 
spacelike region, these take the form 



1 r 



ImF(t')df 
t'-t 



(4) 



with t = 4M 2 (9M 2 ) the pertinent isovector (isoscalar) 
threshold and M„ is the charged pion mass. The basic 
quantity in this relation is the spectral function Im F(t) 
that parameterizes all physical effects that contribute to 
the nucleon form factors. The most general form of the 
spectral function allowed by unitarity is a sum of con- 
tinua and poles. This is exactly the form of the spec- 
tral function we use in our analysis. The low mass con- 
tinua (27r, KK, pm) are included exactly, whereas higher 
mass continua are approximated by effective vector me- 
son poles. 

In our DR approach, the complete isoscalar and isovec- 
tor parts of the Dirac and Pauli form factors, respectively, 
are parameterized as 



aj , a l + b l (l-t/c t )- 2 



m v — t 



pm = E 

V—Vi,V2,-- 



m v -t 2(1 - t/di) 



, (5) 



where i — 1,2. Each pole term comes, in principle, with 
three parameters - the two residua and the mass. While 
the low-mass pole terms can be interpreted as physical 
vector mesons, i.e. the oj and the </>, the higher mass poles 
are effective poles that parameterize unknown continuum 
contributions. The residua for the light vector mesons u), <p 
are related to more or less known coupling constants, see 
e.g. Ref. [7] for a detailed discussion. Therefore, we con- 
strain the residua of the light isoscalar vector mesons as: 
0.5 GeV 2 < K| < lGeV 2 , \a%\ < 0.5 GeV 2 IM and 
\4\ < 2 GeV 2 , \4\ < 1 GeV 2 14 . Additional poles below 
1 GeV 2 can be excluded since no other vector mesons exist 
in this region. At higher masses, the widths get broader 
and with growing distance to the space-like region their 
position cannot be resolved precisely. Therefore this part 
of the spectral function can only be described by effec- 
tive pole terms. Besides the vector meson pole terms, our 
DR approach includes the two-pion continuum, the KK- 
and p7r-continua. The two latter contributions to the form 
factors are also parameterized by pole terms with fixed 
masses and residua, as was done before, see Refs. [15,16, 
I14j . The last term in Eq. ([5]) is an appropriate parameter- 
ization of the two-pion continuum, that also includes the 
p-pole as discussed above. The parameters are: 

oi = 1.084, a 2 = 5.800, h = 0.079, b 2 = 0.751, 
ci = 0.300 GeV 2 , c 2 = 0.225 GeV 2 , 
dt = 0.522 GeV 2 , d 2 = 0.562 GeV 2 . 

Note that this parameterized form of the low-mass con- 
tinua is only used for faster numerical evaluations. The 
total number of isoscalar and isovector poles is determined 
from the stability criterion of Ref. |17j , which in a nutshell 
can be described as using the lowest numbers of parame- 
ters that are required to obtain a good fit to the data. The 
superconvergence relations enforcing the correct asymp- 
totic power behavior of the form factors, 



/ Im Fi(t)t n dt = 0, 

J t 



* = 1,2, 



(G) 



with n = for F\ and n = 0, 1 for F 2 , are included as 
well. 

With this DR ansatz, we reconstruct the differential 
cross sections and fit to the MAMI data for ep scatter- 
ing and simultaneously to the world data for the neutron 
form factors. The fits are done in the one-photon approx- 
imation via the Rosenbluth formula where the data con- 
tains most of the higher order corrections [3]- However, 
the Coulomb corrections are replaced in our analysis by 
a better approximation as was done by |18j . For a recent 
discussion on this point, see Refs. [19114] . We remark that 
two-photon exchange corrections are expected to be small 
for the charge radius extraction, see eg. Refs. [2THI2"!"] . but 
somewhat more sizeable for the magnetic radius [4"ll22|. 
This requires further study. 

We fit to the Mainz data floating the normalization 
of the individual data sets by at most ±4% (in most 
cases this normalization change is less than 1%), using 
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Fig. 1. Form factors from our dispersion relation analysis in 
comparison to world form factor data as given by [£| (updated 
to include also the more recent data [31] V 

the 31 normalization parameters given in Ref. [2 7) . Em- 
ploying the aforementioned stability criterion, we include 
5 isoscalar and 4 isovector resonances with in total 15 ad- 
ditional parameters. We impose the constraints from the 
normalizations, from the superconvergence relations and 
the fact that the masses of the lowest isoscalar poles, the 
w(782) and the 0(1020) are known. For details of the fit- 
ting procedure, see Ref. [7]. This fit gives a good descrip- 
tion of the proton cross section and neutron form factor 
data with x 2 cd = 2.2. The form factors arising from this 
fit are shown in Fig. [TJ The corresponding fit parame- 
ters are given in Tab. [TJ The obtained residua of uj and 
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0.783 


0.747 


0.003 


Vl 


1.205 


0.826 


-2.557 


<i> 


1.019 


0.221 


-0.014 


V2 


2.140 


-0.502 


0.438 


Si 


1.638 


0.977 


0.234 


V3 


1.000 


-0.381 


-0.342 


S2 


2.400 


-0.501 


-0.077 


V4 


1.193 


-0.227 


0.831 


S3 


1.033 


-0.540 


-0.333 











Table 1. Fit parameters from dispersion analysis (cf. Eq. ((5])). 
Masses my are given in GeV and couplings aj in GeV 2 . 



cj) are in agreement with their couplings to the nucleon, 
as described above. All remaining poles have residua with 
natural size. Varying the number of pole terms in both 
isospin channels does not affect the radii significantly but 
is included in our error estimate. 

To address the issue of the theoretical uncertainties of 
our analysis, we vary the isovector and isoscalar contin- 
uum contributions in all possible combinations and repeat 
the fit to obtain the corresponding ranges. The two-pion 
continuum is varied by 5% below the first minimum of the 
pion- nucleon scattering amplitudes and 20% above it. We 
note that at present a Roy-Steiner machinery is set up to 
improve the representation of the t-channel pion-nucleon 
phase shifts required here |23j . but for our investigation 
we use an updated version of the two-pion continuum rep- 
resentation from Ref. [24] , including new data for the pion 
vector form factor from KLOE 25 and BABAR [33] • The 
error in the two-pion continuum from the difference be- 



tween the KLOE and BABAR data sets is well below the 
assigned 5-20% uncertainty. It is also important to note 
that the Roy-Steiner analysis of Ref. [23] has shown that 
the KH80 solution used here is internally consistent. The 
KK- and p7r-continua are varied by 20%. We have also 
varied the number of effective pole terms and obtained 
stable results. 

The extracted proton radii are: 

^ = 0.84™fm, ^=0.861°;" fin. (7) 

For the neutron magnetic radius, we obtain rjf = 0.88 ± 
0.05 fm , consistent with previous analyses. The neutron 
electric radius is very insensitive to the variation of the 
continua. We obtain the value (r^) 2 = —0.127 fm 2 . To 
examine the influence of the neutron data, we also fit the 
DR model exclusively to the electron-proton scattering 
cross sections. The effect on the proton radii is less than 
half a percent on the numbers given above. If the super- 
convergence relations, Eq. ©, are not enforced, the same 
radii are obtained within the quoted uncertainties. 

To assess the stability of the results obtained, we make 
use of a continued fraction (CF) approach as advocated, 
e.g., by Sick [18] . Every rational function can be repre- 
sented by a general continued fraction function and in 
contrast to a polynomial, this can generate poles. There- 
fore it can well approximate the data beyond the threshold 
where the imaginary part sets in. We fit the following in- 
verse CF ansatz [18] . 

F(t) = (8) 



to our data, in order to examine under which circum- 
stances and how exactly the generation of isolated poles 
suffices to simulate the so important two-pion continuum 
contribution. In particular, we have investigated the de- 
pendence of the fits on the highest value of Q 2 (denoted 
Qmax from here on) for which data is included. This de- 
pendence can be used to demonstrate the need to include 
imaginary parts in the fit functions and include data be- 
yond the two-pion threshold to obtain a stable fit. The 
electric radius derived from such fits shows only small 
fluctuations for varying Q 2 lax , as shown in Fig. [5] For 
the magnetic radius, we do not find stable results with 
varying number of continued fractions and thus will not 
consider it any more. Remarkably, the electric radius of 
the proton show large fluctuations in the CF approach 
for Q 2 nax < 0.1 GeV 2 , but then converges quickly towards 
values consistent with the ones based on the theoretically 
preferable DR approach. These results might at first ap- 
pear surprising, as naively a low Q^ax might be favored, 
because the relevant information for the radius extraction 
is related to the region of small momentum transfer. On 
the other hand, the uncertainty due to the normalization 
has obviously a stronger influence on the radius for lower 
Q 2 . The less stable values at very low Qn iax in the DR 
approach might be explained by this normalization un- 
certainty or simply by an insufficient amount of data. At 
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Fig. 2. Dependence of the proton electric radius r v E on Q„iax 
for CF fits. Here, n labels the number of terms. For comparison, 
the DR approach is also shown. 



higher values of Q^ ax the fluctuations in the derived radii 
increase. We assume that this is due to an insufficient 
number of free fit parameters for more data since the fluc- 
tuations are stronger for lower n. Increasing the data range 
also increases the x 2 e d slightly. Of course, this could either 
show the inability of the models to describe the data at 
higher momentum transfers or it could be due to a bad 
estimation of errors in this data. We cannot separate these 
effects, which may both apply here: (i) it is obvious that 
the fit functions are only approximate and (ii) the errors 
of the data are scaled to yield a Xi- C( i ~ 1 for a spline fit 
and thus contain a model-dependence [27] . Therefore, the 
X 2 values or confidence limits are not used to weight the 
extracted radii. The results are weighted equally and only 
the spikes of the fits with x 2 ed — are omitted, less than 
3 for each set of results. The remaining x 2 cd values for all 
fits with Qmax > 0-1 GeV 2 vary, e.g., for the CFs of order 
5 between x 2 ed = 1-03 and 1.64. The average and standard 
deviation of the radius values extracted for Qmax between 
0.1 GeV 2 and the complete data set range up to 0.98 GeV 2 
are r p E — 0.85 ± 0.01 fm for n = 5, n = 8 and n — 12 
(to obtain these results, we have utilized the improved 
Coulomb corrections as described before). These numbers 
are in good agreement with the values obtained via the DR 
approach and the recent muonic hydrogen measurements. 
We are presently studying the application of the confor- 
mal mapping technique of Hill and Paz [9] to the MAMI 
data to further sharpen our conclusion on the value of the 
proton charge radius [3"2"] . 

In this Letter, we have reconsidered the proton radii 
determination from the new MAMI electron-proton scat- 
tering data, using two very different methods. First, we 
have used a dispersive representation of the nucleon form 
factors that includes the so important constraints from 
unitarity and analyticity - these were not considered in 
the fit functions utilized in Refs. [3JH]. Including also data 
from the world data set on the neutron form factors, we 
achieve a good description of the cross section data. We 
have also performed an analysis of the systematic uncer- 



tainties by allowing for generous variations in the 7T7t, KK 
and 7rp contributions to the spectral functions. The re- 
sulting electric and magnetic radii, cf. Eq. ([7]) are in stark 
contrast to the ones obtained in Refs. OH], but our results 
are consistent with all earlier dispersive analyses of the nu- 
cleon form factors [28,29,30,7 . Also, the small value for 
r p E is consistent with the recent muonic hydrogen mea- 
surement pp. We have also applied the continuous frac- 
tion method to the data. By construction, it can emulate 
narrow poles but to account for the so important two- 
pion continuum, one has to choose a sufficiently large fit 
range, Q^ax > 0.1 GeV 2 , to achieve stable results. Again, 
this method leads to values for the proton charge radius 
in agreement with the dispersive approach. We conclude 
that the small proton electric radius is indeed favored if 
analyticity and unitarity are properly included into the 
description of the nucleon electromagnetic form factors. 
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